Loading...
 

Wprowadzenie do rachunku wariacyjnego

Początki rachunku wariacyjnego sięgają drugiej połowy XVII wieku, kiedy to J. Bernoulli sformułował tzw. problem brachistochrony: wyznaczyć tor łaczący dwa punkty tak, aby punkt materialny ześlizgujący się po nim pod wpływem siły ciężkości przebył go w najkrótszym czasie. Sformułujemy go w języku matematycznym.


Załóżmy, że punkt materialny o masie \( \hskip 0.3pc m\hskip 0.3pc \) znajduje się w punkcie \( \hskip 0.3pc A.\hskip 0.3pc \) Mając punkt \( \hskip 0.3pc B\hskip 0.3pc \) leżacy pod linią horyzontalną przechodząca przez punkt \( \hskip 0.3pc A,\hskip 0.3pc \) wyznaczyć linie łączącą \( \hskip 0.3pc A\hskip 0.3pc \) i \( \hskip 0.3pc B\hskip 0.3pc \) tak, aby ześlizgujący się pod wpływem siły ciężkości punkt materialny przebył ją w najkrótszym czasie. Przyjmujemy przy tym, że nie występuje siła tarcia.


Niech \( \hskip 0.3pc A=(0,0),\hskip 0.3pc \) \( \hskip 0.3pc B=(a,b),\hskip 0.3pc \) gdzie \( \hskip 0.3pc a>0,\hskip 0.3pc \) \( \hskip 0.3pc b<0.\hskip 0.3pc \) Zgodnie z prawem zachowania energii mamy

\( \dfrac 12 mv^2- \dfrac 12 mv_0^2= \displaystyle\int_0^ymg ds=mgy,\quad y\in [0,b], \)

gdzie \( \hskip 0.3pc m\hskip 0.3pc \) - oznacza masę punktu, \( \hskip 0.3pc v\hskip 0.3pc \) - prędkość, \( \hskip 0.3pc v_0\hskip 0.3pc \) - prędkość początkową, \( \hskip 0.3pc g\hskip 0.3pc \) - przyspieszenie grawitacyjne. Przypomnijmy, że prawo zachowania energii mówi iż suma energii kinetycznej i energii potencjalnej jest stała. Wzór ( 1 ) wyraża fakt, że zmiana energii kinetycznej jest równa zmianie energii potencjalnej. Z ( 1 ) wynika, że

\( v^2=v_0^2+2gy. \)


Niech \( \hskip 0.3pc x=x(t),\hskip 0.3pc \) \( \hskip 0.3pc y=y(t)\hskip 0.3pc \) będzie szukanym równaniem krzywej. Prędkość punktu wzdłuż tej krzywej określona jest zależnością \( \hskip 0.3pc v(t)= (x^\prime (t),\,y^\prime (t)).\hskip 0.3pc \) Oczywiście możemy przyjąć, że \( \hskip 0.3pc x^\prime(t)>0.\hskip 0.3pc \) Wówczas

\( \vert v(t)\vert=\sqrt{\big(x^\prime(t)\big)^2+ \big(y^\prime(t)\big)^2} =\sqrt{1+\Big(\tfrac{y^\prime(t)}{x^\prime(t)}\Big)^2} x^\prime(t) =\sqrt{1+(y^\prime(x))^2}\dfrac{dx}{dt}, \)

gdzie

\( y(x)=y(t(x)),\hskip 1pc \dfrac{dy}{dx}=\dfrac{dy}{dt}\cdot \dfrac{dt}{dx}= \dfrac{y^\prime(t)}{x^\prime(t)}. \)

Przekształcając ostatnie wyrażenie i uwzględniając równość

\( \vert v(x)\vert=\vert v(t(x))\vert=\sqrt {v_0^2+2gy(t(x))}=\sqrt {v_0^2+2gy(x)} \)

otrzymamy

\( \dfrac {dt}{dx}= \dfrac {\sqrt{1+(y^\prime(x))^2}}{\vert v(x)\vert} = \dfrac{\sqrt{1+(y^\prime(x))^2}}{\sqrt{v_0^2+2gy(x)}}. \)


Całkując ostatnie wyrażenie w przedziale \( \hskip 0.3pc [0,\,a]\hskip 0.3pc \) otrzymamy wzór na czas potrzebny do przebycia krzywej \( \hskip 0.3pc AB\hskip 0.3pc \)

\( T=\displaystyle\int_0^a\dfrac{\sqrt{1+(y^\prime(x))^2}}{\sqrt{v_0^2+2gy(x)}}dx, \qquad \textrm{przy czym}\quad y(0)=0,\hskip 0.3pc y(a)=b. \)

Jak widać czas ten zależy od funkcji \( \hskip 0.3pc y.\hskip 0.3pc \) Problem polega więc na znalezieniu minimum wyrażenia

\( {\rm min}\, T(y) \)

na zbiorze krzywych regularnych \( \hskip 0.3pc y=y(x),\hskip 0.3pc \) łączących punkty \( \hskip 0.3pc A\hskip 0.3pc \) i \( \hskip 0.3pc B.\hskip 0.3pc \)


Po jakiej krzywej powinien ześlizgiwać się punkt materialny startujący z punktu \( \hskip 0.3pc A(0,0)\hskip 0.3pc \) aby w najkrótszym czasie osiągnąć linie \( \hskip 0.3pc x=a.\hskip 0.3pc \)

Z poprzednich rozważań wynika, że należy znaleźć minimum wyrażenia

\( {\min}\displaystyle\int_0^a\dfrac{\sqrt{1+(y^\prime(x))^2}}{\sqrt{v_0^2+2gy(x)}}dx, \)

na zbiorze krzywych regularnych \( \hskip 0.3pc y=y(x),\hskip 0.3pc \) wychodzących z punktu \( \hskip 0.3pc A\hskip 0.3pc \) (tzn. \( \hskip 0.3pc y(0)=0\hskip 0.3pc \)) i przecinających prostą \( \hskip 0.3pc x=a.\hskip 0.3pc \)


Wyznaczyć na płaszczyźnie figurę o największym obszarze mając zadany obwód.


Załóżmy, że brzeg obszaru opisany jest krzywą \( \hskip 0.3pc x=x(t),\hskip 0.3pc \) \( \hskip 0.3pc y=y(t),\hskip 0.3pc \) \( \hskip 0.3pc t\in [\alpha ,\beta ].\hskip 0.3pc \)
Wiadomo, że długość krzywej wyraża się wzorem

\( L=\displaystyle\int_{\alpha }^{\beta }\sqrt{(x^\prime)^2+(y^\prime)^2}dt, \)

a pole powierzchni wzorem

\( S=\dfrac 12 \displaystyle\int_{\alpha }^{\beta }(xy^\prime-x^\prime y)dt. \)

Problem polega więc na znalezieniu maksimum funkcjonału (3) na zbiorze wszystkich regularnych krzywych zamkniętych na \( \hskip 0.3pc [\alpha ,\beta ],\hskip 0.3pc \) spełniających warunek (2).


Powierzchnia o najmniejszym polu przechodząca przez daną krzywą. Niech \( \hskip 0.3pc \Gamma\hskip 0.3pc \) będzie krzywą zmkniętą w \( \hskip 0.3pc \mathbb R^3.\hskip 0.3pc \) Znaleźć powierzchnię \( \hskip 0.3pc S\hskip 0.3pc \) o najmniejszym polu, której brzegiem jest krzywa \( \hskip 0.3pc \Gamma.\hskip 0.3pc \)


Załóżmy że krzywa \( \hskip 0.3pc \Gamma\hskip 0.3pc \) opisana jest równaniami \( \hskip 0.3pc x=x(t),\hskip 0.3pc \) \( \hskip 0.3pc y=y(t),\hskip 0.3pc \) \( \hskip 0.3pc z=z(t),\hskip 0.3pc \) \( \hskip 0.3pc t\in[\alpha ,\beta ].\hskip 0.3pc \) Niech \( \hskip 0.3pc D\hskip 0.3pc \) będzie rzutem szukanej powierzchni na płaszczyznę \( \hskip 0.3pc Oxy,\hskip 0.3pc \) czyli obszarem ograniczonym krzywą \( \hskip 0.3pc x=x(t),\hskip 0.3pc \) \( \hskip 0.3pc y=y(t),\hskip 0.3pc \) \( \hskip 0.3pc t\in[\alpha ,\beta ].\hskip 0.3pc \) Należy znaleźć funkcje \( \hskip 0.3pc z=u(x,y),\hskip 0.3pc \) \( \hskip 0.3pc (x,y)\in D,\hskip 0.3pc \) realizującą minimum funkcjonału

\( \displaystyle\iint_D\sqrt{1+u_x^2+u_y^2}\,dxdy \)

i spełniającą warunek \( \hskip 0.3pc u(x(t),y(t))=z(t)\hskip 0.3pc \) dla \( \hskip 0.3pc t\in[\alpha ,\beta ] .\hskip 0.3pc \)


Załóżmy, że dana jest powierzchnia regularna \( \hskip 0.3pc S\hskip 0.3pc \) oraz punkty \( \hskip 0.3pc A\hskip 0.3pc \) i \( \hskip 0.3pc B\hskip 0.3pc \) leżące na tej powierzchni. Znaleźć krzywą o najmniejszej długości łączącą punkty \( \hskip 0.3pc A\hskip 0.3pc \) i \( \hskip 0.3pc B,\hskip 0.3pc \) całkowicie leżącą na powierzchni \( \hskip 0.3pc S\hskip 0.3pc \) (tzw. krzywą geodezyjną).


Przyjmijmy, że powierzchnia \( \hskip 0.3pc S\hskip 0.3pc \) zadana jest równaniem

\( r=r(u,v)=\big(x(u,v),y(u,v),z(u,v)\big),\qquad (u,v)\in D. \)

Ponadto przyjmijmy, że szukana krzywa ma równanie

\( r=r(u(t),v(t)), \qquad t\in[\alpha ,\beta ]. \)

Z analizy wiadomo, że długość krzywej (4) wyraża się wzorem

\( \displaystyle\int_{\alpha }^{\beta } \sqrt {E(u^\prime)^2+2Fu^\prime v^\prime+G(v^\prime)^2}dt, \)

gdzie

\( E=r_u\cdot r_u,\hskip 0.3pcF =r_u\cdot r_v,\hskip 0.3pc G= r_v\cdot r_v, \)

(symbol \( \hskip 0.3pc \cdot \hskip 0.3pc \) oznacza iloczyn skalarny, \( \hskip 0.3pc r_u=(x_u,y_u,z_u), \hskip 0.3pc \) \( \hskip 0.3pc r_v=(x_v,y_v,z_v).\hskip 0.3pc \))

Problem polega zatem na znalezieniu minimum funkcjonału (5) na zbiorze krzywych regularnych

\( u=u(t),\hskip 0.4pc v=v(t),\hskip 0.5pc t\in[\alpha ,\beta ], \)

takich że: \( \hskip 0.3pc \big(u(t),v(t)\big)\in D\hskip 0.3pc \) dla \( \hskip 0.3pc t\in[\alpha ,\beta ],\hskip 0.3pc \) \( \hskip 0.3pc r\big( u(\alpha ), v(\alpha )\big)=A,\hskip 0.3pc \) \( \hskip 0.3pc r\big(u(\beta ), v(\beta )\big)=B.\hskip 0.3pc \)


Ostatnio zmieniona Środa 06 z Grudzień, 2017 18:42:15 UTC Autor: Julian Janus, Józef Myjak
Zaloguj się/Zarejestruj w OPEN AGH e-podręczniki
Czy masz już hasło?

Hasło powinno mieć przynajmniej 8 znaków, litery i cyfry oraz co najmniej jeden znak specjalny.

Przypominanie hasła

Wprowadź swój adres e-mail, abyśmy mogli przesłać Ci informację o nowym haśle.
Dziękujemy za rejestrację!
Na wskazany w rejestracji adres został wysłany e-mail z linkiem aktywacyjnym.
Wprowadzone hasło/login są błędne.